Microdroplet impact at very high velocity 

Claas Willem Visser Q Yoshiyuki Tagawa, Chao Sun[j] and Detlef Lohs^] 

Physics of Fluids Group, Faculty of Science and Technology, 
J. M. Burgers Centre for Fluid Dynamics, 
University of Twente, The Netherlands 
(Dated: March 7, 2013) 

Abstract 

Water microdroplet impact at velocities up to 100 m/s for droplet diameters from 12 to 100 \xm 
is studied. This parameter range covers the transition from capillary-limited to viscosity-limited 
spreading of the impacting droplet. Splashing is absent for all measurements; the droplets always 
gently spread over the surface. The maximum spreading radius is compared to several existing 
models. The model by Pasandideh-Fard et al [lj agrees well with the measured data, indicating 
the importance of a thin boundary layer just above the surface, in which most of the viscous 
dissipation in the spreading droplet takes place. As explained by the initial air layer under the 
impacting droplet, a contact angle of 180 degrees is used as model input. 
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I. INTRODUCTION 



High-pressure spray cleaning, droplet-wall interactions in diesel engines, and plasma 
spraying are notable examples of processes in which high-speed impact of small droplets 
on a solid surface is a key phenomenon. In these applications, droplets with a characteristic 
size of 1 to 100 |xm and a velocity of order 100 m/s impact on a solid surface [2]-[4]. Despite 
this industrial interest, microscale droplet impact at very high velocity {Uq > 50 m/s) has 
only been studied for solidifying metal droplets [5]. This is mainly due to the challenging 
parameter regime: very high spatial and temporal resolutions are required to study the rel- 
evant phenomena. In addition, it is difficult to create impact events at these velocities. An 
understanding of the phenomena of high-speed micro-sized droplet impact is thus lacking 

In this work, we aim to extend current results in three ways. First, we present a 
novel high-velocity droplet generation, by using a method to create ultrafast liquid jets. 
Second, using high-speed imaging, the impact dynamics is studied. Third, a quantitative 
investigation of the maximum spreading radius will be presented and compared to existing 
models, to improve our understanding of droplet spreading. 



II. PARAMETER SPACE 

To compare our results to previous work, a phase diagram of the droplet size (Do) and 
impact speed (Uq) is plotted in figure [j] (a). Most work up to now has focused on impact 
dynamics of droplets with a size of ^1 mm in diameter. The studies on microdroplet impact 
were mainly at relatively low speed (up to 10 m/s)[5j EJ Ej. Our study connects these 
previous investigations, in particular those of water microdroplet s at lower velocies [8] and 
metal microdroplet impact at very high velocities [5]. 

Figure [j] (b) shows a phase diagram of the achieved Reynolds and Weber numbers of 
experimental droplet studies. The Weber number is defined as We = pD^U^/a, where a 
the surface tension and p the density. The Reynolds number is given by Re = pD Uo/p 
where p represents dynamic viscosity. As shown in figure [j] (b) , the Re- We values of our 
measurements largely overlap with previous data. However, our data have been taken for 
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microdroplets instead of mm-sized droplets as indicated in figure [j] (a). As the droplet size 
is a key control parameter for impact dynamics, various models developed to describe the 
mm-sized droplet impact may not hold for microdroplet impact dynamics. The large impact 
velocity also explains why our data have relatively high Weber numbers for given Reynolds 
numbers. 

Our data cover the transition between the capillary- to the viscosity-dominated limits 
of droplet spreading [TO] , as shown by the solid line in figure [j] (a) (an explanation of 
this transition is provided below). So far, to investigate this transition, liquids of different 
viscosities and surface tension were required in order to achieve a sufficient coverage of the 
Re- We parameter space. Using droplets of microscopic scales lowers the Weber number for 
which this transition takes place, allowing to study the transition region with a single liquid. 

III. EXPERIMENTAL SETUP 

To create high- velocity microdroplets (U Q > lOm/s), we make use of a new method to 
create ultrafast liquid jets [18], as sketched in figure [2] (a). In a nutshell, the method works 
as follows. By focusing a laser pulse with a microscope objective, a vapor bubble is created 
in a capillary tube, by laser-induced cavitation. From this bubble, a shock wave travels to 
the meniscus, which in turn forms a liquid jet thanks to flow focusing. Subsequently, this 
jet breaks up into tiny droplets with a velocity similar to the jet velocity (figure [5] (b)). 

It is found that the tube diameter is a key control parameter for both the jet diameter 
and the jet tip velocity [18]. Therefore, capillaries with diameters are used to create a range 
of droplet sizes and velocities. In addition, the laser energy and the distance between the 
laser focus and the meniscus are varied to generate droplets at different velocities for a given 
tube diameter [18j, resulting in the diameter- and velocity ranges, as shown in figure [2] (c). 
With this method, approximately 170 droplet impacts on a dry surface were examined. For 
the sake of clarity, the data is binned in the figure. 

To create droplets at velocities between 1 and 10 m/s, a commercially available Microdrop 
dispenser is used. By varying the input voltage between 60 V and 160 V per pulse, a range 
of velocities is covered. The droplet diameter ranges from 40 \xm to 80 |Lim, as shown in 
figure [2] (c). 

A standard microscope slide was used as an impact plate, which is placed above the tip of 
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FIG. 1: Parameter space of (a) droplet velocities and radii and (b) Weber and Reynolds numbers 
for various droplet impact experiments [T} 15} IHHT7] . The solid blue line signals P — 1 and separates 
the capillary regime with P < 1 from the viscous regime with P > 1. 



the capillary tube. Atomic force microscope measurements indicated a roughness R a below 
10 nm. The droplet impact is visualized from the side. After each impact measurement, 
the glass plate was cleaned with ethanol and dried with paper tissue. Frequent checks with 
a 10 x optical microscope (after the cleaning procedure) indicated that this method usually 
resulted in an optically clean surface, i.e., hardly any paper fibers or dirt were sticking to the 
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FIG. 2: (a) Setup used for the generation of fast droplets, (b) Jet generation and breakup at 
different instants, (c) Range of the achieved droplet velocities and radii, for different droplet 
generator settings. Capillary tube diameters are decreased from 500 fim to 200 (im to generate 
smaller and faster droplets. In addition, to create slow droplets, a MicroDrop apparatus is used. 
As approximately 230 measurements were performed, the data are binned. The bars represent one 
standard deviation. 

surface. The liquid used in the present work was a standard water-based blue inkjet printer 
ink, with density p = 998 kg/m 3 , surface tension a = 72 mN/m 2 , and viscosity /i = 10 3 
Pa-s. 

Magnifications of 10 x to 40 x are obtained by combining a standard Olympus 10 x 
objective with an adjustable 12 x zoom lens (Navitar 1-50015). An Olympus ILP-1 light 
source is used for illumination. An ultra high-speed camera (Shimadzu HPV-1) is used 
to study the microdroplet impact dynamics, at recording rates of 1.25 • 10 5 to 10 6 frames 
per second. For the fastest droplets (Uq « 100 m/s), the impact duration is approximately 
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r = Dq/Uo ~ 0.2 |lls, which is below the temporal resolution of the camera. Thus, com- 
pletely capturing of such events requires even higher frame rates [19] or pulsed illumination, 
e.g. as used in Ref. [20]. However, these techniques require a level of control of the moment 
of impact which is not achieved with the present setup. The high-speed camera has a 
minimum shutter time of 500 ns. Therefore, at high impact velocities, substantial motion 
blur is observed (e.g. figure [3] (c)), as the shutter time approaches r. Still, with the current 
setup, the presence of splashing could be assessed and the maximum spreading radius could 
be observed in an entirely new parameter regime. 



IV. RESULTS 

A. Impact Phenomenology 

In general, the droplet impact process can be divided into the following phases. When 
the droplet approaches the solid surface, the air between the falling drop and the surface is 
strongly squeezed, leading to a pressure buildup in the air under the drop. The enhanced 
pressure results in a dimple formation in the droplet and an air layer development between 
the droplet and the target plate [2TH23 J. Before the droplet wets the surface, the liquid 
moves on top of this air cushion. The droplet extends in the radial direction until it reaches 
a maximum spreading radius. In this phase, splashing can occur [71 [T6l [2T1 [24H26] . Finally, 
the droplet completely wets the surface and an air bubble is entrapped [23j . 

The first aim of this work is to assess the dynamics of microdroplet impact, i.e. to find out 
whether the drop is in the splashing or gentle spreading regime. The latter one is defined by 
droplet deformation into a pancake shape, without satellite droplet formation. Time series 
of droplet impacts and the subsequent spreading phases are shown in figure [3j For these 
figures, the velocity ranges from 0.7 m/s to 100 m/s. 

At low velocity (Uq = 0.7 m/s), initial flattening of the droplet bottom is observed, as 
shown in figure [3] (a). Subsequently, the droplet spreads over the substrate into a (virtually) 
half-dome shaped cap. During spreading, the droplet starts to oscillate and comes to rest 
only after the droplet has reached its maximum diameter. Additionally, air bubble entrain- 
ment is observed (the small black dot, just left of the droplet center). These phenomena are 
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FIG. 3: Time series of droplet impact at different impact velocities, (a) Uo = 0.7m/s, Dq = 70(am, 
We = 0.52, Re = 56, and the time lag St between the frames is 8 (as. Droplet oscillations and 
air bubble entrapment can clearly be seen, (b) Uq = 7.7m/s, Dq = 71 (am, We = 60, Re = 613, 
and St = 4 (as. Spreading into a thin center sheet and a thicker rim are observed, followed by 
oscillations and partial withdrawal of fluid from the rim towards the droplet center, (c) Uo = 73m/s, 
Do = 23(am, We = 1.8 • 10 3 , Re = 1.9 • 10 3 , and St = l(as. The details of the spreading phase can 
no longer be resolved, (d) Uq = lOOm/s, Do = 20(am, We = 3 • 10 3 , Re = 2.3 • 10 3 , and St = l\is. 

consistent with what had been reported in Ref. [8], which includes a detailed discussion on 
the droplet oscillation frequency and the size and cause of the bubble/cushioning entrained 
[T6l [221 [231 [27] . At medium velocities (7.7 m/s, figure [3] (b)), the droplet deforms into a 
disc-like structure. Here, the central impact area is a sheet-like structure surrounded by a 
thicker rim. Again, capillary oscillations were observed. Finally, even at very high velocities 



7 



(figures [3] (c) and[3](d)), still no splashing is observed. Thus, we conclude that gentle impact 
occurs for all velocities and droplet sizes investigated in the present work. 

B. Maximum spreading 

Now the maximum spreading radius will be determined and compared against various 
models. As knowledge of the maximum spreading diameter is of paramount importance 
for industrial applications, a plethora of models has been developed [U [9j [TUl EHH32J. It 
is generally agreed on that the spreading is limited by either viscosity or surface tension. 
Therefore, a key issue is to define which of these is dominant. To study this issue, we will 
first briefly summarize several models. 

Assuming an inviscid liquid, the maximum spreading is limited by the surface tension. 
Balancing the Laplace pressure force with the inertial deceleration of the drop, a scaling 
of D max /D ~ We 1/4 is obtained P HO]. This scaling is remarkably robust [TQl fT4l [T6] in 
the capillary regime. As shown in figure [i] (a), the data of mm-sized droplets impact on 
superhydrophobic surfaces obtained by Tsai et al [16J well agree with this 1/4 scaling law. 

Another limiting case is to completely ignore surface tension and assume that the max- 
imum spreading radius is limited by the viscous dissipation during droplet spreading [TO] . 
This yields D max /D$ ~ Re 1 / 5 as scaling law, which holds well for mm-sized droplets in the 
viscous regime as shown in Ref. [10J. 

To quantify the transition between the viscous and the capillary regimes, Clanet et al 
[TO] defined the parameter P = We/Re 4 / 5 . For P < 1, a surface-tension dominated regime is 
expected (i.e. D max /D$ ~ We 1 / 4 ), whereas a viscous scaling is predicted for P > 1, yielding 
the previously mentioned D max /D ~ Re 1 / 5 . In figure^ (a), the transitional Weber number 
We t) defined by P — 1, for droplets of D = 2 mm and 50 |xm are plotted. For 2 mm-sized 
droplets, We t ~ 3 x 10 3 , whereas it decreases to We t ~ 2 x 10 2 for 50 |Lim-sized ones. As 
shown in figure [4] (a), the present microdroplet data are in between the transition regime 
from the capillary regime to the viscous regime. For a given We, the microdroplet spreading 
is lower than that of mm-sized droplets due to the viscous effects. This combination of 
two simple models provides a decent first description of the impact dynamics. However, 
more detailed models are possible and have been developed. As described by Chandra & 
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Avedisian [29j, the dissipated energy equals the work done, and can thus be estimated by 



U 



2 



W = j j 0«Kr0«Kr/x^J (1) 

where (j) the dissipation function, estimated as /ll(Uo/L) 2 jTUl [29] , with r the typical impact 
timescale, V v = TrLD^ax ^ e total droplet volume, and L the characteristic dissipation 
length scale, which in this model is selected as the height of the droplet splat h. In the 
model by Chandra & Avedisian [29J, W is used in an energy balance Ek + E s — W + E' s , 
where Ek is kinetic energy, and E s and E' s are surface energies before and after impact, 
respectively. However, as shown by Pasanideh-Fard et al [1] and our results, this model 
strongly overpredicts the maximum spreading radius. 

The model was revisited by Pasandideh-Fard et al [lj. Their numerical simulations 
suggested that the dissipation is taking place in a thin boundary layer within the expanding 
droplet, implying that the characteristic length scale L in the above eq. [T]is not the pancake 
thickness h but has to be replaced by the Prandtl-Blasius boundary layer thickness 5 [I]: 

(2) 



'Re 

finally resulting in the following equation for D max /D \ 



D max /D = J ^ (3) 

' y 3(1 - cos 0) + 4(We/ V^) 

Here, 9 is the contact angle. At low impact velocities (We < 10) this model saturates, as 
the droplet impact can be considered to be effectively static and Ek and W vanish. 

In figure [4] (b) , the measured maximum spreading for the microdroplet impacts is plotted 
versus the Weber number. The maximum spreading is measured at the moment when the 
deformation of the droplet is maximum before it wets the surfaces, as shown in figure [3] 
(a,b). At low velocity, the present data saturates around a spreading of ^1.3 times the 
initial diameter, which has a good overlap with previous low-speed microdroplet impact 
experiments [8j. At high velocities, spreading of only ^3 times the initial diameter is found 
for even the fastest droplets. This is consistent with (even faster) impact of molten metal 
microdroplets [5], even though the liquid properties were very different from our experiments. 

Figure [i] (b) is complemented with the models of Chandra & Avedisian [29j and 
Pasandideh-Fard PQ, evaluated for an initial diameter of D = 50|xm. We first examine 
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FIG. 4: The normalized maximum spreading D max /Do vs. We. (a) Solid circles: the present 
microdroplet data. Open squares: mm-sized droplets impact on superhydrophobic surfaces from 
Tsai et al [16j . Dash-dot line: the capillary model by Clanet et al jTU] • Solid line: the 
viscous model by Clanet et al. [10J. The short vertical dashed lines show the transitional Weber 
number We t at P — 1 for a 2 mm-sized droplet and for a 50 \im droplet, (b) Solid circles: the 
present microdroplet data. Stars: data from Mc Donald et al [5] with impact of molten metal 
microdroplets. Open triangles: low-speed microdroplet impact experiments by Van Dam & Le 
Clerc [8J. Lines: Model by Chandra & Avedisian [29] (dark-green lines) and by Pasandideh-Fard 
P] model (black lines), evaluated for the initial diameter Dq = 50(am and contact angles tv/Q 
(dashed lines) and tt (solid lines) 



the maximum spreading calculated from these models for two different contact angles, i.e. 
6 = 7r/6 and n. The dashed lines in figure [i] (b) clearly show that the contact angle of tt/6 
results in a much larger spreading factor for both models as compared to the experimental 
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data. As discussed above, an air layer is present between the spreading drop and the solid 
surface before the droplet wets the surface [2TH24] , implying a "contact angle" which would 
be best described by 6 = tt. As shown with the solid lines in figure [3] (b), this approach re- 
markably decreases the deviations between the experiments and the results of both models. 
A good agreement is found between the microdroplet data and the Pasandideh-Fard model 
[lj up to We ~ 10 3 . This indicates the importance of a finite boundary layer thickness in 
the dissipation of spreading droplets in the present parameter regime. For We > 10 3 , the 
increasing trend of D max /D versus We seems to saturate. This finding is consistent with 
Ref. [9j, and will be studied in future work. 

V. CONCLUSIONS 

The impact of water microdroplets on a smooth solid surfaces is investigated experimen- 
tally. By using a new droplet-generating device, impact events were created at velocities 
from 1 to 100 m/s and droplet diameters between 12 and 100|Lim. This parameter regime 
covers the transition between surface tension- and viscosity-dominated spreading of the 
droplet. For all impact events, no splashing is observed. The maximum spreading radius 
was compared to several models. The model by Pasandideh-Fard et al [1\ performs best, 
indicating that boundary layer dynamics play a key role in droplet spreading. In addition, 
we find that an initial contact angle of 180 degrees should be used as input value. This 
confirms the presence of an air layer under the impacting droplet. 
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